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Abstract 

The cascade mass reconstruction approach was used for mass re- 
construction of the lightest stau produced at the LHC in the cascade 
decay g — > bb — ► ^bb — > f\rbb — > Xi T1 ~bb. The stau mass was re- 
constructed assuming that masses of gluino, bottom squark and two 
lightest neutralinos were reconstructed in advance. 

SUSY data sample sets for the SU3 model point containing 160k 
events each were generated which corresponded to an integrated lu- 
minosity of about 8fb _1 at 14 TeV. These events were passed through 
the AcerDET detector simulator, which parametrized the response of 
a generic LHC detector. The mass of the f\ was reconstructed with a 
precision of about 20% on average. 



* Permanent address: Institute for Nuclear Research, 60-th Oct. pr. 7a, Moscow 117312, 
Russia 
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I. Introduction 



If supersymmetry exists at an energy scale of 1 TeV, SUSY particles such as 
gluinos and squarks should be abundantly produced at the LHC. Assuming 
R-parity conservation, these particles cascade down to the lightest supersym- 
metric particles (LSPs). A detailed discussion of possible SUSY effects at the 
LHC is given in [I]. 

In this paper we consider a mass reconstruction of the lightest stau (fi) 
in the cascade decay 




Figure 1: A gluino cascade decay chain. 



The study of sleptons and squarks of third generation is of a special 
interest. Their masses can be very different than that of sparticles of the 
first and second generation, because of the effects of large Yukawa and soft 
couplings in the renormalization group equations. Furthermore they can 
show large mixing in pairs (t^, (b^, 6r) and (f^, tr). The properties of f\ 
are also important for determination of the dark matter relic density. 

The reconstruction of a SUSY event is complicated because of escaping 
LSPs and many complex and competing decay modes. At present there 
are two different approaches to SUSY mass reconstruction. The "endpoint 
method", which has been widely studied [2] - [9] for the LHC at high in- 
tegrated luminosity of about 100 — 300 fb _1 , looks for kinematic endpoints 
of invariant mass distributions. The second method of SUSY particle mass 
reconstruction is the "mass relation approach" [10] - [UJ , based on the mass 
relation equation which relates SUSY particle masses and measured momenta 
of detected particles. It was shown in [12] that the "mass relation approach" 
can be used for luminosities as low few fb . 

In this work the cascade mass reconstruction approach [12] , that we de- 
veloped earlier and applied to a different SUSY cascade decay, is used for 
mass reconstruction of the lightest stau at the LHC energy of 14 TeV with 
an integrated luminosity of about 8 fb -1 . The stau mass is reconstructed 
assuming that the masses of the gluino, bottom squark and two lightest neu- 
tralinos were reconstructed in advance. At high integrated luminosity the 
stau mass could be extracted directly from the endpoint of rr invariant mass 
distribution. However at integrated luminosities below 10 fb -1 the problem 
becomes quite challenging because there is no a sharp edge in the rr invariant 
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mass distribution because of escaping neutrinos, a reconstructed stau mass is 
very sensitive to uncertainties in input parameters (neutralino masses) and 
there is a high level of SUSY background. 

The cascade mass reconstruction approach [1 2] applied for stau mass 
reconstruction in this work is based on a consecutive use of the endpoint 
method, an event filter and a combinatorial mass reconstruction method. 
The endpoint method is used to get a rough estimate of the stau mass and 
the corresponding errors. This first estimate of mass is used by an event 
filter for each event in the data sample to reduce background for the follow- 
ing mass reconstruction step. Finally, the stau mass is reconstructed by a 
maximization of a combined likelihood function, which depends on all five 
sparticle masses (gluino, bottom squark, stau and two lightest neutralinos) , 
and is constructed for each possible combination of five events in the data 
sample. 

II. Simulation 

We choose for this study the SU3 model point. This point has a significant 
production cross section for the chain (JTJ; gluinos and squarks should be pro- 
duced abundantly at the LHC The bulk point SU3 is the official benchmark 
point of the ATLAS collaboration and it is in agreement with the recent 
precision WMAP data [13]. This model point is described by the set of 
mSUGRA parameters given in Table (CQ). 



Point 


m 


mi/2 


A 


tanfi 


fj, 


SU3 


100 GeV 


300 GeV 


-300 GeV 


6 


> 



Table 1: mSUGRA parameters for the SU3 point. 



Assumed theoretical masses of SUSY particles in the cascade (JTJ), the 
total branching ratio and a cross section generated by IS A JET 7.74 [TJJ are 
given in Table (T2]). 



Point 




m b\ 


m h 2 


m x° 


m tau 1 




BR 


a[pb] 


SU3 


720.16 


605.93 


642.00 


223.27 


151.46 


118.83 


2.71% 


19 



Table 2: Assumed theoretical masses of sparticles, branching ratio BR and 
production cross section a at the SU3 point. Masses are given in GeV. 



Branching ratios for the gluino decay chain (jJJ) at the SU3 point are 



~ 9.2% 7 16.6% ~n 4 8-7% ~ 100% ~n n ^r(W 

g — > b 2 — > x 2 — * r i — ► Xi °- 75 % 

Monte Carlo simulations of SUSY production at model points were 
performed by the HERWIG 6.510 event generator [T5]. The produced 
events were passed through the AcerDET detector simulation [TH], which 
parametrized the response of a detector (LHC detector descriptions can be 
found in [IT], |18j). The efficiency for jet reconstruction and labeling was 
80%. An additional factor of 50% took into account tau efficiency, which 
typically corresponded to a rejection factor of about 100 for QCD jets back- 
ground. Samples of 160k SUSY events were used. This approximately corre- 
sponds to 8 fb _1 of integrated luminosity for the SUSY SU3 point production 
cross section of 19 pb at 14 TeV. Five different sets of 160k SUSY events were 
considered to demonstrate the stability and precision of the mass reconstruc- 
tion approach. 

As input parameters for the stau mass reconstruction procedure, hypo- 
thetical masses of the gluino, bottom squark and two lightest neutralinos are 
used. These input masses are assumed to have been reconstructed by the 
cascade mass reconstruction approach for SUSY data sets corresponding to 
4 fb _1 P2] and are presented in Table El 



Set 


m- g 


m~ h 


m x° 




1 


701±57 


600±57 


208±21 


98±22 


2 


712±55 


608±53 


254±21 


143±20 


3 


664±78 


564±80 


219±24 


109±23 


4 


767±62 


649±65 


258±35 


148±34 


5 


655±45 


545±47 


208±21 


96±20 



Table 3: Reconstructed SUSY particle masses and reconstruction errors for 
five data sample sets. 

In order to isolate the chain ([1]) the following cuts were applied: 

• two jets tagged as r with opposite charge satisfying transverse momen- 
tum cuts pt > 30 GeV and > 25 GeV 

• two b-tagged jets; 

• at least four jets, satisfying p T1 > 100 GeV, px2 > 50 GeV, px3 > 
50 GeV, p T4 > 50 GeV; 

• M eff > 500 GeV and EJp iss > 0.2M e// , where EJp iss is the event's 
missing transverse energy and M e f f is the scalar sum of the missing transverse 
energy and the transverse momenta of the four hardest jets; 

• invariant mass of r-tagged jets satisfying 10 GeV < M TT < 300 GeV. 
Note that at the first stage of the reconstruction procedure, the stau mass 

is estimated by considering the chain —> X2Q ~* TT 2 q — > Xi TT( l with the 
same cuts as mentioned above except requiring b-tagged jets. 
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It was shown in [6] that the Standard Model processes are suppressed sig- 
nificantly by the above requirements on r-tagged jets. The Standard Model 
dominant backgrounds surviving the hard cuts are Z + jets and ti production, 
where both W's decay leptonically into a bbll state. This background was 
estimated to be 1/10 of the SUSY backgrounds [6]. The ti background can 
be reduced to about 2% by applying the event filter procedure of the cascade 
mass reconstruction approach [12J. Thus, the remaining Standard Model 
contribution is negligible in comparison with that of SUSY backgrounds, 
and therefore Standard Model background is not included in the following 
analysis, the SUSY background is about three times the signal rate. 

Table (0J shows the number of signal events and SUSY background events 
for the SU3 model point after cuts were applied to the five sets of 160k SUSY 
events. The dominant SUSY background consists of r + and r~ produced in 
two different decay chains. The classification of events as signal and SUSY 
background is based on simulated truth information. The SUSY background 
to the process ([I]) is significant. It follows from Table (j3J) that for the SU3 
point the number of SUSY background events is a factor 3 greater than the 
number of signal events. 



Set 


Total 


Signal 


SUSY Backg. 


Ratio 


1 


199 


47 


152 


4.2 


2 


212 


56 


156 


3.8 


3 


216 


53 


163 


4.1 


4 


216 


48 


168 


4.5 


5 


229 


65 


164 


3.5 


1-5 


1072 


269 


803 


4.0 



Table 4: The number of signal and SUSY background events after cuts ap- 
plied to 160k SUSY events. Ratio = (Signal+Background) /Signal. Row 1-5 
shows the number of events for five sets combined. 



III. Preliminary estimate of stau mass by the 
endpoint method with low statistics 

The endpoint method has been widely used to determine masses of SUSY 
particles [3], in particular it has been applied to the decay chain [6], [T9j - 

QL -> X°<? hrq -> Xi^rq (2) 

which is a subprocess of the cascade (jTJ if one considers q\ instead of b. 

Fig. (T2]) shows the rr invariant mass distribution for Set 1 of 160k events 
at the SU3 point after application of the kinematic cuts and after subtraction 
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of same sign r-tagged (background) jets was performed. This allowed to 
reduce background because in background processes positive and negative 
jets are mostly uncorrelated. 
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Figure 2: Invariant mass distribution of r-tagged jets for Set 1. 

The distribution in Fig. (j2j) does not have a sharp edge because of un- 
detected neutrino energy in tau-decays but the end-point is clearly seen at 
about 109 GeV and can be reconstructed. 

To extract endpoints from the rr invariant mass distribution a fit to this 
distribution was performed using a parabola plus a straight line that matched 
the parabola. 

The results of the fit for the five data sets are shown in Table (jHJ). 





Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


endpoint 


109.3±1.1 


105.0±1.0 


104.0±0.7 


108.5±1.0 


106.8±0.6 



Table 5: Endpoints determined from fitting of rr-edges. 



These results have to be compared with the theoretical endpoint of 101.7 
GeV for the SU3 point. The fit was performed by MINUIT [22J which gave 
errors presented in Table (jSJ). These uncertainties are small because of a 
specific form of the fit function. More conservative estimates of uncertainties 
[6], [20] are about 4 GeV but in our case this contribution is not important 
because the most significant uncertainties come from 20% uncertainties in 
sparticle masses. 

Once endpoints are found the stau mass can be preliminary estimated 
from the analytical expression for the rr endpoint 

m 2 TT =(i- f)(f-x)/f (3) 
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The following notations are used for masses y = m? , f = m% , £ = m 2 . Q . 

The results of the fit for preliminary stau mass estimates and their errors 
based on found endpoints, edge reconstruction errors, sparticle masses and 
their errors are summarized in Table ([6]). 



fx 


Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


solution 1 


149±135 


209±48 


174±56 


205±151 


159±51 


solution 2 


137±137 


173±47 


138±54 


187±150 


125±50 



Table 6: Stau masses preliminary determined from fitting of rr-edges. For 
each set two stau mass solutions are given. 

It is seen from Table (JHJ) that uncertainties of stau mass reconstruction 
are very large and additional steps in stau mass reconstruction are required. 

Eq. gives two solutions for stau mass. Because we use these masses as 
a rough estimate for more precise procedures described below in the following 
an average of these two solutions is used as an estimate for the stau mass. 
Final results are not sensitive to this approximation. In particular, a good 
estimate for the stau mass at this stage would be just an average of two 
lightest neutralino masses. 

IV. Background suppression 

As a second step of stau mass reconstruction an event filter is used to suppress 
background before the final fit. At the event filter stage we assume that three 
light SUSY particle masses (x®, fx, Xx) are fixed- The mass of fx is taken as 
the preliminary mass found by the endpoint method in the previous chapter. 
Neutralino masses are taken from Table [3J Gluino and sbottom masses are 
assumed to be distributed uniformly in the range rrig ± 2a, m~ b ± 2a where 
masses and errors are given in Table [31 

The event filter procedure is based on minimization for each event of the 
function 

X 2 K, m b ) = — + A/(rM (4) 

i=x Gi 

where index i runs over two b-quarks and two taus and labels the mea- 
sured absolute momenta p™ eas , the uncertainties <7j in their measurement, 
and the event true absolute momenta ^" ent . The approximate function ([41 
takes into account only uncertainties in tau and b-jet energy measurements. 
Note that positions of each of two b-jets and of each of two taus in the decay 
chain ([T]) are unknown. It is quite simple to resolve the b-jets assignment 
because usually (about 96% of the time) the b-jet with higher p T originates 
from the 6-quark decay. Therefore, for each event we assume that the b-jet 
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with higher px originates from the 6-quark decay. Taus with higher px (after 
cuts) are produced in both vertices with comparable probabilities. In this 
work it is assumed that taus with higher pt originate from \2 decay. We use 
the following parametrization for <7j in equation (Tj0): for b-jets and tau-jets 
a/E = O.b/yfEjGeV) © 0.03. The Lagrange multiplier A in Eq.® takes 
into account the mass relation constraint f(rh,p) [10] - [12] which relates 
masses and momenta of particles in the chain ([1]) . 

For signal events, the event likelihood distribution has a maximum in 
the region of the (g, b) mass plane correlated with the true masses of g 
and b. Thus signal events should give a peak in the region of true masses. 
For background events there is no strong correlation of maximum likelihood 
distribution with true (g, b) masses. Therefore if we chose arbitrary a point in 
the (g, b) mass plane in the range m g ± 2<r, m~ b ± 2a x 2 reaches its minimum 
in this region with much higher probability for a signal event than for a 
background one. For each event, 10 5 points are generated randomly in the 
mass plane range rrig ± 2a, m- b ± 2a and the \ 2 is calculated. If x 2 < 10 in 
at least 1000 points this event is considered SIS db SI gnal candidate and it is 
retained for the subsequent analysis. 

After the application of the event filter the ratio of background events 
to signal events is reduced approximately by a factor 1.5 as can be seen 
from Table [7J The contribution of the combinatorial background is given ap- 
proximately by ((Signal+Background)/Signal) 5 , corresponding to five event 
combinations required at the last step of mass reconstruction, and is therefore 
significantly suppressed. The suppression factor varies from 5 to 10. Note 
that any five event combination including at least one background event is 
considered background. 



Set 


Total 


Signal 


SUSY 


Ratio 


Suppression 




Number 


Events 


Backg. 




Factor 


1 


199/98 


47/33 


152/65 


4.2/3.0 


5 


2 


212/118 


56/49 


156/69 


3.8/2.4 


10 


3 


216/116 


53/39 


163/77 


4.1/3.0 


5 


4 


216/114 


48/37 


168/77 


4.5/3.1 


6 


5 


229/124 


65/55 


164/69 


3.5/2.3 


8 


1-5 


1072/570 


269/213 


803/357 


4.0/2.7 


7 



Table 7: The number of signal and background events before/after an appli- 
cation of event filter to 160k SUSY events. The last row shows the sum over 
all five sets. Ratio = (Signal+Background) /Signal. The last column gives 
approximately a suppression factor (Ratio 5 before to Ratio 5 after) for five 
event background combinations. 
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V. Stau mass reconstruction 



As a third step a combinatorial procedure [32] is used for the final stau mass 
reconstruction. It is applied only to the events that pass the event filter. 
At the final stage of mass reconstruction when the physical background has 
already been reduced, we will consider all possible five event combinations 
from the event sample. 

SUSY particle masses are reconstructed by a search for a maximum of 
a combined likelihood function constructed for each possible combination of 
five events in the data sample. 

The x 2 function for an event is defined by 



(movent r ,meas\2 jj I ^ event „_ \2 

XL* = £ - 7* ' +E ' " 2 " + A,/ + A 2 /" . (5) 

' * (J. ^— ' (T 

i=l 1 n=l n 

where the first term takes into account deviations of measured momenta of 
b-jets and tau-jets from the true ones. The second term takes into account 
that the masses of sparticles vary from event to event as approximated by 
a Gaussian of width a n instead a Breit-Wigner distribution. In Eq.(J5]) the 
mass relation and rr edge constraints are taking into account by Lagrange 
multipliers Ai, A2. Standard deviations corresponding to the mass widths are 
taking to be 15 GeV for the gluino, 5 GeV for bottom squark and 1 GeV for 
light masses. The first two numbers are comparable with theoretical widths 
for heavy SUSY particles. The last number takes into account the fact that 
light SUSY particles are quite narrow or stable. We note that the results 
of the mass reconstruction are not strongly sensitive to the actual values of 
sparticle widths. 

At this step it is assumed that the starting distributions of the heavy 
masses are uniform in the range: mean value ± 2a, where the mean values 
and standard deviations are given in Table [3J For the three light masses a 
Gaussian distribution is assumed with mean values and standard deviations 
found by the endpoint technique as given in Table [6j The MINUIT code is 
used to search for the minimum of the x 2 function for five event set with five 
mass parameters. 

The reconstructed SUSY particle mass distribution and a result of fitting 
this distribution by a Gaussian are presented in Figj3]for Set 1. As can be 
seen in this figure, the reconstructed mass distribution is well described by a 
Gaussian except for small tails. The mean values of the Gaussian fit are the 
reconstructed sparticle masses. 

Final results of this mass reconstruction approach are presented for five 
data sample sets of 160k events each in Table [HJ 

It follows from Table [8] that reconstruction errors after the final step are 
significantly reduced in comparison with preliminary estimates [3] obtained by 
end-point method. 
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Figure 3: A reconstructed stau mass distribution for Set 1, including SUSY 
background with integrated luminosity 8 fb _1 . The line is a result of Gaussian 
fit. 



Particle 


Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


fx 


154±24 


192±23 


168±25 


202±29 


150±22 



Table 8: Reconstructed stau masses and reconstruction errors for five data 
sample sets of 160k events each. 



In order to illustrate a spread in reconstructed masses the results of Table 
[8] are also shown in a form of ideogram [23] in FigJD for the five data sample 
sets. Each reconstructed mass in an ideogram is represented by a Gaussian 
with a central value error crj and area proportional to 1/Vj. The solid 
curve is a sum of these Gaussians. 

The Gaussian form of the ideogram and relatively small shift of peak 
positions with respect to theoretical masses demonstrate the self-consistency 
of stau mass reconstruction approach. 

IX. Conclusion 

We applied the cascade mass reconstruction approach developed in [T2] to 
reconstruction of f\ mass at the LHC with a low integrated luminosity of 
about 8 fb at 14 TeV. This luminosity can be reached in the early stage 
of LHC operation. At low integrated luminosity stau mass reconstruction is 
complicated because escaping neutrinos modify the edge of the rr invariant 
mass distribution, reconstructed stau mass is very sensitive to uncertainties 
in input parameters (neutralino masses) and the level of SUSY background 
is high. 
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Figure 4: An ideogram of reconstructed stau masses, including SUSY back- 
ground for five data sample sets with integrated luminosity 8 fb _1 . The 
triangle marker gives the position of the SU3 theoretical mass. Points with 
error bars correspond to five data sample sets. 

Our approach to the stau mass reconstruction was based on a consecu- 
tive use of the endpoint method, an event filter and a combinatorial mass 
reconstruction method. 

The endpoint method allows preliminary estimate of f\ mass by assuming 
that sparticle masses required as input parameters for this procedure are 
known from the same experiment with a precision of about 10% for heavy 
sparticles (such as gluino and bottom squark) and about 20% for the two 
lightest neutralinos. In this work these input masses were obtained by the 
cascade mass reconstruction approach [12J. At high integrated luminosity of 
100 - 300 ftr 1 they could be found by using one of techniques proposed in 

The stau mass was reconstructed with a precision of about 20% 
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